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In this paper we review a recent proposal to understand the long time limit of glassy dynamics 
in terms of an appropriate Markov Chain. [1]. The advantages of the resulting construction are 
many. The first one is that it gives a quasi equilibrium description on how glassy systems explore 
the phase space in the slow relaxation part of their dynamics. The second one is that it gives an 
alternative way to obtain dynamical equations starting from a dynamical rule that is static in spirit. 
This provides a way to overcome the difficulties encountered in the short time part of the dynamics 
where current conservation must be enforced. We study this approach in detail in a prototypical 
mean field disordered spin system, namely the p-spin spherical model, showing how we can obtain the 
well known equations that describes its dynamics. Then we apply the same approach to structural 
glasses. We first derive a set of dynamical Ornstein-Zernike equations which are very general in 
nature. Finally we consider two possible closure schemes for them, namely the Hypernetted Chain 
approximation of liquid theory and a closure of the BBGKY hierarchy that has been recently 
introduced by G. Szamel. From both approaches we finally find a set of dynamical Mode-Coupling 
like equations that are supposed to describe the system in the long time/slow dynamics regime. 


I. INTRODUCTION 

Glass is an out of equilibrium state of matter. Let by themselves, the microscopic configurations and macroscopic 
properties of glasses slowly change in a process known as physical aging, where the free-energy evolves towards lower 
and lower values. This process is very slow and is characterized by a separation of time scales, where fast degrees 
of freedom appear to be in thermal equilibrium in the background of the slow degrees of freedom. The evolution of 
the latter is sometime described as an effective random walk in configuration space, depicted as a rough free-energy 
landscape, where the system wanders from one metastable state to another. The rules by which the metastable 
states are selected in the dynamical process determine the observed properties of the systems. A leading role in the 
comprehension of slow off-equilibrium dynamics of glassy systems is played by mean field theory (see [2] for a review), 
which allows to describe asymptotic aging regimes through scaling laws and effective temperatures associated to the 
violations of the fluctuation dissipation theorem [3, 4]. The analysis performed in [5, 6] and in [7-11] related these 
effective temperature to quasi-equilibrium selection of metastable states during glassy dynamics. This notion of quasi¬ 
equilibrium exploration can be formalized through the introduction of a suitable Markov Chain where fast times are 
effectively coarse grained, and it is assumed the set of states available at any give time, which are the ones at a specific 
distance from the present state, are selected by a Boltzmann law [1]. This chain construction reproduces the results 
of long time relaxational dynamics within mean-field theory and it gives the right long time dynamical equations 
whenever the system has a finite configurational entropy. We believe that it captures the principles of exploration of 
configuration space in glasses also in realistic systems where metastable states are sufficiently long-lived. Moreover, 
since the short timescales are completely coarse grained, the method automatically produces time reparametrization 
invariant equations. Thanks to the Boltzmann prescription, the chain share many formal features with equilibrium 
systems. This observation straightforwardly suggests simple ways to treat long time dynamics of structural glassy 
systems taking advantages of standard approximation schemes originally devised to study equilibrium [12]. Moreover, 
since short times are coarse-grained, constraints such as energy and mass conservation that complicate short time 
analysis are not relevant here. In this contribution we review the main results of our approach and present some new 
results and derivations that were just hinted in previous publications. The paper is organized as follows: in the first 
section we introduce the basic dynamical construction and we discuss the response properties of the system and its 
equilibrium measure. In the second section we discuss spin glass mean field theory, and we give a new derivation of the 
dynamical equations through a probabilistic analysis. This avoids the use of replicas which were used in a previous 
publication [I]. Then we present the results of the direct integration of our equations that were never published before. 
In the third section we review the applications to realistic liquid models. We present the derivation of a dynamical 
version of the Ornstein-Zernike equation. We complement this equation with two closure schemes, that give both a 
final dynamical equation that is of the same kind of standard Mode-Gouling Theory (MCT). Within this scheme we 
are able to predict the properties of the dynamics both in the equilibrium regime and in the aging one. 
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II. THE BOLTZMANN PSEUDODYNAMICS CONSTRUCTION 

In this section we review the Boltzmann Pseudodynamics construction recently introduced in [1] and we show some 
generalities about correlation and response functions that can be computed using this formalism. Let us consider a 
system described by a set of internal degrees of freedom that we call Si and that will be addressed as spin variables 
(our notation is very close to the one encountered in spin systems but can also be used to treat particles in a liquid 
where the internal degrees of freedom are the position of the particles). In the following we will define a dynamical 
rule to evolve such system so that we will indicate with Si{t) the configuration of the system at time t. The Boltzmann 
Pseudodynamics (BPD) is a discrete time dynamics defined from the following dynamical rule: given the configuration 
of the spins at time t, the configuration at time t + 1 occur with a probability that is given by 

M{S{t + l)|5(t)) = {cit, t + l)-q {S{t),S{t + 1))) 

Z[0t+i;m]= E + + 

S(t+1) 

The function q(a, r) is an overlap function that measures the similarity between the configurations a and r. For spin 
system it is given simply by q{a, r) = ^i^i/N where Ci and Ti are the values of the spin in the two configurations. 
For particle systems instead one can choose between many different definitions. A simple possibility is to use the same 
function as for spin systems through a lattice gas coding of the liquid configurations e.g. by dividing the volume in 
small cells and defining binary variables that code for occupation numbers of the cells^. The probability of a trajectory 
given an initial configuration A(0) at time t = 0 is given by 


t-i 

P[S{t),S{t - 1),.. .,^(1)15(0)] = P[Sm [] MiSik + 1)|^(A:)) (3) 

fc=i 

where P is a given measure over the inital conditions. Note that in the above dynamics we fix the set of variables 
{C{t + I,t)} and {/?(} from the outside. In glassy dynamics the temperatures are naturally fixed to the value of 
the thermal bath, while C{t + l,t) should be chosen in a self-consistent way in order to achieve the right separation 
between fast and slow time scales. While in certain applications it can be interesting to consider a time dependence 
of the temperature, from now on in this paper we will consider the case of /? constant in time. We will show in the 
following that this pseudodynamics provides a coarse grained description of real time dynamics in which fast processes 
are seen as instantaneous. For finite values of the choosen values C{t + 1, t), at each time step then the system chooses 
a configuration at a macroscopic distance from the previous one. This configuration is chosen as an equilbrium one 
at the prescribed distance. In this sense the “fast” time scales that in real dynamics are needed to equilibrate within 
a metastable state are coarse-grained. The main assumption is that all configurations satisfying the constraint q = C 
are equally reachable by the fast relaxation processes. This assumption is fine if C is properly chosen: e.g. a value 
close to the typical overlap qsA- 

A computer implementation of (1) would require for each step of the chain a Montecarlo simulation in which the 
fast time scale is reinserted to achieve the prescribed sampling. For that reason we call it pseudodynamics. We will 
be interested to the long chain limit in which the relevant time scales are much larger than the unit of the elementary 
time step and the system moves to even larger distances than the ones reached in a single step. 


A. Response functions 

A fundamental characterization of a glassy dynamics is provided by the linear response functions. Their relation 
with fluctuations during aging shows, both in mean-field theory [3, 13] and in simulations [14] the emergence of 


^ Another popular choice in systems of particles is to use 

= ^^w{\xi-yj\) ( 2 ) 

id 

where X = {fci, ...,Xn} and Y = {yi, ...,yAr}, are two configurations of N particles and where w is some positive short ranged function, 
e.g. iu(r) = exp(— 
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effective temperatures ruling the exchanges of energy between slow degrees of freedom. It has been proposed in [15] 
that non trivial response and effective temperatures in non-equilibrium dynamics are possible if trajectories (clones) 
starting in the same point and subject to different thermal noise get separated during their evolution. One can define 
the clone correlation function of a simple observable as for example the magnetization m{S^(t)) as follows 

Qs{t,u) = Es'(s) [E(TO(t)|5(s))E(TO(u)|^(s))] for t,u> s. 

The internal expectations denoted by E(-) are the averages over independent thermal trajectories that start from the 
same initial configuration ^{s) at time s, while E 5 ( 5 ) is the average over these initial configurations. In [15] it was put 
forward the conjecture that non trivial effective temperatures, are only possible for systems such that for large time 
separation Qs(t,u) tends to the same minimal value as the correlation C{t,u) itself. This is different from systems 
undergoing domain coarseing in phase separation where Q remains much larger then C and a non trivial response in 
the aging regime is absent. Unfortunately despite empirical evidence, a formal relation between response and clone 
correlation function was lacking. In the dynamics we just introduced this relation emerges naturally for pseudo-times 
M = s -I- 1 and t > u corresponding to real times t ^ u and s,u such that C{u,s) is in the beta relaxation regime. 
Consider the dynamics (1) in a time dependent field ht coupled with an observable m{S(t)), function of the system 
configuration ^(t). The Hamiltonian in presence of the field is 

= (4) 


The response function is defined as usual 


A{t,s) 


d{m{S{t))) 

dhs 


( 5 ) 


where the average is done over the multiple realizations of the trajectories of the system. Because of the causal 
structure of the Markov Chain (I), the response function is non zero only if t > s. To analyse this quantity we start 
from 

- 1), = P im{S{s)) - E[m(5(s))l^(s -l)])P - 1), ...,^(I)|^(0)) (6) 

where 

E [m(^(s))15(s - I)] = z[i3;s\s-l)] ^ {q{Sl,S{s - 1)) - C{s, s-l)). (7) 

This leads to 


A(t,s) = P[{'m{S{t))'m{S{s))) - {m{S{t))E[m{S{s))\S{s - 1)))] = P[C{t,s) - D{t,s)\ (8) 

where 

(7i(^(s))) = E • • • E Ms{s))p -1), ^(1) \m) Hm), 

m s(o) 

C{t, s) = {m{S{t))m{S{s ))), D{t, s) = {m{t)E{m{S'{s))\S(s - 1))) = Qs-i^t, s). (9) 

Notice that, by simple properties of conditional probability, if t < s one has D{t,s) = C{t,s) and the response is 
A(t, s) = 0 as it should be. Eq. (9) provides the announced relation between response and clone correlation: a 
response at large times is only possible if Qs-i{t, s) = D{t, s) differs from C{t, s). 

We would like now to make some remarks on the long chain (i.e. long times) limit. By properly choosing the 
constraining value close to the typical overlap value, C « Qea, all the fast dynamics is coarse-grained in a single step of 
the pseudodynamics. Consequently, we have C{t, t) « qea and the equal time response A(t, t) = — (m(t))^], 

corresponding to the equilibrium-like response in the beta regime, is expected to be a quantity of order one, which 
coincides with /3[C'(t, t) — C{t + l,t)] since D{t,t) should be close to C{t + l,t). On the other hand if t > s (to 
be intended as t s in real times) one can expect D{t,s) to be close to C{t,s). For a chain of length t the total 
susceptibility to a field that acts from time I to t, is x(I) = A{t, s). This quantity remains hnite for t —)• oo if 
A(t, s) = R{t, s)ds is infinitesimal in the continuous time limit. In this way the sum converges to 

X{t) = A{t,t) + [ dsR{t,s), (10) 

Jo 

where the first and second terms are respectively the responses from the fast and the slow dynamics. A non trivial 
response in the long time regime is thus associated with a non zero response function R{t, s) i.e. to decaying clone 
correlation function that in the continuum limit is Qs{t, s + ds) = C{t, s) — dsTR{t, s), since in the fast dynamics the 
fluctuation-dissipation relation dsC{t, s) = T R{t, s) holds. 
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B. Equilibrium measure 

Despite we are interested to the applications of (1) to glassy dynamics and time scales where equilibration does not 
occur, in general, for time independent correlations (^(s + l,s) = C and finite system’s volumes, the Markov chain 
(1) is ergodic and it is interesting to study its equilibrium measure. This is not the ordinary Boltzmann distribution. 
In fact we can observe that the detailed balance is verified with respect to the modified distribution 

= ( 11 ) 

^2 


where 


Za = ^ (q{S, ff) - C) (12) 

S,S' 

This is therefore the equilibrium distribution of the chain. Equivalently one can see that the measure for two config¬ 
urations at consecutive times is 


= Ae-/5[^(S)+^^(s')]j (q{S,S') - C) . 


(13) 


III. MEAN-FIELD GLASSY DYNAMICS 


In this section we would like to analyse mean field spin glasses and show how to obtain a full characterization of the 
dynamics in terms of a single correlation function and its conjugated response function. The analysis was previously 
performed in [1] with the aid of replica method that can be employed to treat analytically the denominators in 
(1). Here we propose an alternative derivation that avoids the use of replicas similar to the derivation of mean-field 
dynamical equations for Langevin dynamics presented in [16]. 

We consider specifically the spherical p-spin model which provides the canonical example of mean-field glassy 
dynamics. The Hamiltonian of the model is 

N N 

Hj[S;h] = - = P[J,,„.iJcxexp 

ii<...<ip i=l i=l 

where we have introduced a site dependent magnetic field in the system that is needed in order to compute the 
response function. In order to make the presentation as simple as possible we will restrict our analysis to the case 
p = 3 even if the general case is a straightforward generalization of this one. 

Due to the mean field nature of the model, we can get closed dynamical equations in terms of two point correlation 
and response functions. Let us consider an arbitrary function or operator (/)(S) dependent on a trajectory S = 
{S[(u)}u=o s-iid write the obvious identity 

E./ J [<^{S)P{S\S{m = 0 (15) 

where, Ej represents the average over the disordered couplings. In order to obtain equations for correlation and 
response functions it is enough to consider the insertion of 4>{S) = Silr) and (j){S) = which do not depend on 

the quenched variables J, we can therefore average direcly the measure P(S|5(0)) 


NP- 


-Jt 


p! 


(14) 




Ej 


1 

2 


P'^J^jkSj{a)Sk{cr) + Ha S^{a) 


+,^aS,{a - I) + i^a+i[S,{<J + 1) - E {S,{a + l)|S(a))]] P(S|5(0))} . 


(16) 


In order to simplify our analysis we suppose that ^(0) is chosen randomly with uniform probability on the sphere 
and P(S[(0)) is independent of J. In the case where P{S_{0)) depends of J additional terms would appear (see e.g. 
[16]). We can now integrate by part the Jijk in the first term of (16) which results in the substitution Jijk —>■ 
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F? Ee/3[‘5'*(e)'S'i(e)S'fc(e) -E(5'*(e)S'j(e)S'fc(e)|S(e- 1))] to get 


E, 


dSi{a) 


P{S\m)= Ej 


e j,k 

+fiaS,{a) + v^S,{a - 1) + v.+i[SM + 1) - E (5,(a + l)|S(a))]] P{S\S{0))} ■ (17) 


We will make at this point the crucial hypothesis that for typical trajectories jfJ2j ^'^'1 

ly Ei'S'i('^)E(S'i(cr)|S((T — 1)) are self-averaging quantities and coincide respectively with C{T,a) and I?(t, cr). This 
hypothesis is equivalent to the factorization, for i ^ j, {Si{T)Si{a)Sj{T)Sj{a)) = {Si{T)Si{a)){Sj{T)Sj{a)) and anal¬ 
ogous formulas for more then two indices. This “clustering conditions” of the correlation functions, which imply that 
the averages are dominated by a single pure state, exclude a-priori replica symmetry breaking (RSB) effects. It is of 
course possible to include RSB effects even if in the long chain limit, RSB effects are indistinguishable from violations 
of FDT within the RS formalism. Using the factorization hypothesis we finally get 


E, 


dSi{a) 


P(S|5(0))= Ej 


^/3^/3[5,(e)C(e,a)2-E(^,(e)|S(e-l))P(e,a)2] 


+fi,S,{a) + - 1) + u^+i[Sda + 1) - E (5,(a + l)|S(a))]] P(S|^(0))} . (18) 


Inserting at this point 4’{S) = Si{T) and <('(5') = and performing the sum over the trajectories we get respectively: 


-5(r - cr) = ^ /3 [C'(t, e) C{(t, ef - D{t, e) D{a, ef] P a) + v„C{t, a - 1) + j^^+iA(t, cr -f 1) , (19) 

(7 

0 = ^ /3A(e, r) lC(a, - D{a, e)^] -f cr) + J^<tA(t, a - 1) + 5r,^+iVa+iI\{a + l,a Pl) . (20) 

e—T 

These equations have a causal structure that is inherited from the chain construction and can be integrated iteratively 
step by step. In the next section we show that assuming the existence of a long chain limit, these equations reduce to 
the long time equations for the slow part of the Langevin dynamics of the same model. The consistency of this limit 
can be checked through explicit integration of the equations. At low temperatures, one finds an aging regime where 
selfconsistently C{t, t -I- 1) —)■ qea while —> 0 at large r. 


A. The long chain limit 

The equations (19,20) can be easily generalized to arbitrary spherical long range spin glass models with Hamiltonian 
of p-spin mixture type = Ep®p7Ip[5] and correlation function E{H{S)H{S')) = Nf{S ■ S'/N), with f{q) = 

7 Ep It is interesting to write them in the general case: 

r 

-S{r -cr) = ^/?[C'(T,e) f{C{a,e) - D(r,e) f{D{a,e))] + fj.^C{T,a) + iy^C{r,a - 1) + z/^+iA(T,a + 1) (21) 

e^O 

(7 

0 = '^l3A{e,T) [f'{C{a,e)) - f{D{a,e))] P HaA{T,a) P v^A{T,a - 1) -f 5r,a+iVa+iA{a Pl,aP 1(22) 


In the continuum time limit, using D{t, s) = C(t, s) — T R{t, s)ds for t > s, and fixing vt = 0, that corresponds to 
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using the saddle-point value qsA for C{t,t — 1), we get 


n{t)C{t,u) = 

u) = 
fi{t) = 


p[ ds f'{Cit,s))R{u,s) +p f ds f"{C{t,s))Rit,s)Ciu,s) 
Jo Jo 

+/3"(/'(l) - nqEA))Cit,u)+l3^f{C{t,u))0 - qEA) 

13 j dsf"{C{t,s))R{t,s)R{s,u) 

J U 

+[3f'{C{t, u))R{t, u)il - qEA) + PifO) - nqEA))R{t, u), 

T + (3^{f{l)-f{qEA)) 

+/3 f ds if{C{t,s))R{t,s) + f{C{t,s))R{t,s)C{t,s)) 

Jo 


(23) 


where we used the condition C{t, t) = 1 and C{t, t) = qEA- These are well know equations that in the dynamic theory 
of mean field spin glasses, depending on the models, they have a dynamical phase transition below which there are 
aging solutions where fluctuation dissipation theorem and time translation invariance do not hold. 

In order to check the consistency of the aging solutions and the long time limit of the original Markov Chain, 
we have integrated explicitly the discrete BPD equations for the pure p-spin model for p = 3 with f{q) = ^q^ and 
for a mixture of p = 2 and p = 4 with f{q) = It is interesting to consider both models since, as well 

known, the former has a one-step replica symmetry breaking (IRSB) phase and displays aging with a single (inverse) 
effective temperature [3] (3eff = PX while the latter has a full replica symmetry breaking (fRSB) phase and has a 
continuum set of effective temperatures Peff{q) = PX{q) = f"'{q)/f"{q)^^‘^ that during aging depend continuously 
on the value of the correlation function C{t,u) = q. We integrated the equations in the low temperature regime with 
1) fixed to the theoretical value of qEA found from self-consistency of the long time equations with vt = 0. An 
example of the results obtained from the numerical integration of the equations is given in Fig. 1 where we show the 
correlation function of the p-spin model for p = 3 at a low temperature. A popular way of analyzing aging systems 
and detecting effective temperatures consists in looking at the parametric plot of the integrated response function 
Tx(t, s) = ^('^5 e) —)■ T f* du R{t, u), as a function of the correlation function C{t, s). According to mean-field 

theory, for large times these curves tend to a master curve whose slope should be precisely —X{C). 

In Fig.s 2 and 3 we present such plot respectively for the p = 3 model and for the 2-1-4 model, and compare it with 
the asymptotic result predicted by the aging mean field theory. We see that in both the IRSB and fRSB cases, the 
asymptotic aging limit is approached after a very limited number of steps in the chain. We confirm in this way that 
BPD gives a faithful representation of aging dynamics. 


IV. REPLICAS 

In the previous section we used a path-integral representation of the probability of a trajectory to get a set of 
closed equations for the correlation and response functions. In the case of systems with long range interactions an 
alternative way to produce the dynamical equations is by treating the denominators that appear in the definition 
(1) of the Boltzmann Markov Chain with the replica method. This consists in substituting in all the theoretical 
computations the denominator appearing in (1) by a positive power of the constrained partition function to get 

m + l)l^(t)) = ZIP; (C(t, t + 1) - q (A(t), ^(t + 1))) (24) 

The new chain depends on the parameters rit+i that for each time t counts the “number of replicas”. As usual 
in the replica method, these number are considered integers in intermediate computations, but sent to zero in the 
replica limit where (24) coincides with (I). Renaming the “master replica” S{t) = and introducing rit — l “slave 
replicas” S°'{t), a = 1, ...,nt — I, we get: 

-1 "t+i-i ^ 

Mn,+AS°it+l)\S\t))= e-^^"=o'” Yl s(c{t,t + l)-q{S%t),S‘^{t + l))'^ (25) 

The dynamical equations for the correlations and response functions in the case of the p-spin spherical model (23) 
were first derived with this formalism in [1]. Within the replica method it is natural to introduce the total partition 
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C[t, u] 



FIG. 1: The correlation function C{t,u) in the spherical p-spin model for p = 3 from the explicit numerical integration 
of the equations as a function of u/t < 1 for t = 9{blue),19{orange),29{green),39{red). The integration is performed at 
T = 0.517 <Td = -^/S/S = 0.612. C(t + 1, t) is kept fix at the value qea ~ 0.621. The data seem to indicate that C{t, u) tends 
to an exclusive function of the ratio u/t for large u and t. 

Chi T 



FIG. 2: The plot of rescaled susceptibility Tx{t, u) as a function of C{t, u) in the p-spin model for p = 3. The lines are computed 
at times t = 9{blue), 19{orange), 29{green), 39{red) and the temperature is T = 0.517 {Td = \/3/9> = 0.612). We also plot the 
FDT line 1 — C and the modified FD prediction 1 — gsA + X{gEA — C). The value of the fluctuation-dissipation ratio predicted 
by the long time dynamics of Langevin equation is X = 0.610 and coincides well with the data. 


function up to time t 


Ztot{t)= Y. ( 26 ) 

S“(ti) u,a 

where in the sum and the product u runs up to t and, for each u, a runs from 0 to — 1. This expression, complicated 
as it may be, has the same formal structure of a partition function of a replicated equilibrium system. The difference 







FIG. 3: The plot of the rescaled susceptibility Tx{t,u) as a function of C{t,u) for t = 10{blue),20{orange),30{green) in the 
2+4 model with f{q) = + aq'^) and a = 0.1 at temperature T — 0.34. The critical temperature of the model is Tc = 1 . 

We also plot the FDT line 1 — C and the modified FD prediction from RSB, which show the good agreement of the BPD with 
the expected results for infinite time prediction. 


with the usual case comes in the fact that here we have an explicit interaction between replicas at subsequent times 
due to the chain constraint. As in the usual case however, the finite rit system can be interpreted as a mixture of 
interacting particles of different kind. In this way this is a good starting point for approximations. In fact one can 
apply quite straightforwardly all existing approximations for equilibrium mixtures modulo parametrizations of the 
quantities of interest that allow for the analytic continuation needed to take the replica limit. As usual, considerations 
of symmetry under permutations of different replicas play an important role. In the present case, differently from 
the equilibrium case where all replicas are equivalent, the partition function (26) is only invariant under independent 
permutations of groups of slave replicas with the same time index. 

In replica calculations it is quite simple to see that a prominent role in manipulating (26) is played by the replica 
correlation function Qab{s,u) = {S°'{s)S^{u)). This codes for both the correlation and the response functions of the 
previous section. Assuming replica symmetry (i.e. invariance under the previously mentioned group of permutations), 
we can write the most generic replica symmetric parametrization for the overlap matrix 

Qab{s, u) = C(s, u) + A(s, u) + da.iA(M, s) + (5s,„da.6A(M, u). (27) 

In the long chain continuous time limit this becomes: 

Qab{s, u) = C{s, u) + R{s, u)du 0{s — u)SbX + R{u, s)ds 9{s — u)6a,l + [C{u, u) — C{u, M)](5s,Mda.6- (28) 

In [1] it has been shown that this parametrization reproduces the correct dynamical equations for this model in the 
slow dynamical regime. In the next section we will deeply use this form in the context of the replicated liquid theory 
in order to obtain a set of new dynamical equations for structural glasses [12]. 


V. BOLTZMANN PSEUDODYNAMICS FOR SUPERCOOLED LIQUIDS 

In this section we want to develop the BPD formalism to describe the slow regime of the dynamics of supercooled 
liquids undergoing to a glass transition. As we emphasized in the previous section this can be done quite easily. Indeed 
the BPD construction can be directly applied to the replicated liquid theory that has been successfully employed to 
describe the glass transition [17-22]. In particular we can take the equations of the replicated liquid theory used to 
obtain the statics of structural glasses and plug the pseudo dynamics ansatz inside them. Practically, in doing this 
we have to promote simple replica indices to BPD replica indices a (t,a). 

Let us see how all this procedure works. We start from the definition of the basic quantities that can be treated in 
the theory of the replicated liquid [18]. The simplest objects we need are the density field and its two point correlation 
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function that are defined by 


N 

Pa{x) = (^5{x-x^t'’)) Pap{x]y) = (Y^5 {x - x^“'’)5{jj - xf'^)) (29) 

i=l [ij] 


where the sum over [ij] runs on all i, j if a ^ /3 and over i ^ j li a = j3. Moreover we define 


hap{x,y) 


Pafi{x,y) 

Pa{.x)pp{y) 


(30) 


in what follows we will always look for a uniform solution for the density field such that Pa{x) = p. 

We introduce also the direct correlation function Cap{x) that is defined by the replicated Orstein-Zernike (OZ) 
equation: 


'^tot p 

Cap{x) = hap{x) haj{x - y)pcp^{y). 

7 = 1 “' 

It is convenient also to rewrite the same equations in Fourier space 

ntot 

Capiq) = hap{q) - ^ ha.y{q)pcpj{q). 

7=1 

In the glass phase and in within a IRSB picture, the static version of the OZ equation can be written as 


Hq) = c{q) + p [Hq)c{q) + {m- l)h{q)c{q) 

Hq) = c{q) + p \hiq)c{q) + c{q)h{q) + (to - 2)h{q)c{q) 


(31) 


(32) 


(33) 

(34) 


where to is the number of replicas. In the limit to —>■ 1 we get the OZ equations that are needed to compute the 
dynamical MOT transition point 


h{q) = c{q) + ph{q)c{q) 

Hq) = c{q) + p \h{q)c{q) + c{q)h{q) - h{q)c{q) 


(35) 

(36) 


where c and c are respectively the diagonal and off-diagonal parts of the matrix Cap and the same is true for hap. 
These quantities can be thought as the corresponding correlations functions for the supercooled liquid for times that 
are such that any dynamical correlation function is close to its plateau value. In particular they can be used to 
compute the MCT non ergodicity parameter. 

We want to show that once the relevant generalization of (28) is considered for Cap and hap, the equation (32) has 
formally the structure of a long time Mode Coupling equation which can be used both to describe equilibrium slowing 
down as the glass transition is approached and the aging dynamics below the glass transition point. 

As in the simple case of (28) where there is no space dependence, the replica dependence of h and c encode for 
correlation and response functions and 


hapix) = hab{s, U, x) = h{s, u; x) -I- 6abSsu^h{s, s; x) + SaiRhiu, s; x)ds -I- SbiRhis, u; x)du 
Cap{x) = Cabis, u, x) = c(s, u] x) + Sab^su^c(s, s; x) + SaiRc(u, s] a;)ds -I- SblRcis, u; x)du. 

Plugging these forms inside Eq. (32), we get a dynamical version of the OZ equations 

h{q; s, u) = c{q; s,u) + p [h{q; s, 0)c{q; 0, u) + h{q; s, u)Ac{q; u, u) + Ah{q; s, s)c(q; s, u)-|- 

+ ^ J dzh{q]s,z)Rc{q;u,z) + ^ J dz Rh{q] s, z)c{q; z,u) 

Ah{q; s, s) = Ac{q; s, s) + pAh{q; s, s)Ac(q; s, s) 


Rhiq,s,u) = Rciq;s,u) + p 


1 /'® 

Rh{q;s,u)Ac{q;s,s) + Ah{q]u,u)Rc{q]s,u) +-js / dz Rh{q-, z,u)Rc{q; s, z) 

P J U 


(37) 

(38) 


(39) 

(40) 

(41) 

(42) 
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These equations are not closed: we need to provide some kind of closure scheme in order to have a complete set of 
dynamical equations. This will be done in the following. For the moment, let us investigate the properties of these 
relations. It is quite simple to see that these equations are compatible with time translation invariance (TTI) and 
fluctuation dissipation theorem (FDT): 


h{q; s, u) = h{q; s — u) 

d 

Rh{<l]s,u) = e{s - u)j3 — h(q]s - u), (43) 

and analogous relations for c and Re- If TTI and FDT are inserted in the OZ equations, they reduce to the following 
equation (we can set u = 0 due to TTI): 


h{q] s) = c{q;s) +p 


h{q;s)Aco{q) + Ahoiq)c{q;s) + - / dzRh{q; s - z)c{q; z) + h{q] s)c{q;0) 

P Jo 


(44) 


where we have introduced the following notation 

Ah{q-, s, s) = Ahoiq) = 0) - h{q; 0) Ac{q; s, s) = Aco(g) = c{q; 0) - c{q; 0) (45) 


where h{q; 0) and h{q; 0) are nothing but a solution of the static OZ equations. 


A. Closure schemes 

The OZ relations alone are not enough to write a self-consistent system of equations and we need to provide a closure 
scheme for them. Here we discuss two closure schemes: the first one is the standard Hypernetted Chain Approximation 
(HNC) [23] that has been developed extensively to study glasses [17, 19]. The second one is a closure scheme that has 
been introduced by G. Szamel [24] in order to derive the standard MCT equations for the non ergodicity parameter 
from the replica approach. Both the two approaches have advantages and disadvantages. On the one hand the 
HNC approximation is known to not provide a quantitative sensitive description of the glass transition [22]. On the 
other it is variational since it can be derived from a partial resummation of the diagrams that give the free energy. 
This makes it quite suitable for a systematic improvements on top of it. The Szamel’s closure instead is an ad hoc 
scheme to obtain quantitatively the same non-ergodicity factor of standard MCT from the replica approach. The 
disadvantage is that this procednre is not variational and it needs the external input of the static structure factor as 
it is usual in MCT. However quite remarkably, using the BPD construction we are able to derive from these purely 
static approximation/closure schemes a set of dynamical equations that are nothing bnt MCT equations in the long 
time regime. 


1. Dynamical HNC equation 


The HNC closure equation for a replicated system is given by 

ln[/ia/3 (a^, 2/) + 1] + Idfpap {x,y) = hap {x,y)- Cap (x, y) (46) 


where for us /papixiy) = 4‘aT-,ba{x,y) = ST,aSa,b4>ix — y) + St,^iST,a+iw{x — y)v{a) contains the inter particle potential 
at eqnal time and replica indexes multiplier constraining the value of the overlap at consecutive times. Plugging the 
parametrization (37-38) into (46) we obtain for s ^ u 


\a[h{x] s, m) -I- 1] = h{x; s, u) — c{x; s, u) 

h{x; s, u) 


Rc(x;s,u) = Rh(x;s,u) 


h(x; s,u) + 1 


(47) 

(48) 


The dynamical equations (39-42,47,48) provide a complete set of equations that can be solved in time [12]. 

It is quite evident from the BPD construction that the equations that we can get from it must be covariant under 
time reparametrization. Technically this means that if we have a solution h{q;s,u) and c{q;s,u) for these equations 
we can obtain another solution from this one in the following way: we consider a monotonically increasing function 
f{t); we can write a new solution as 

h'{q; s, u) = h{q; f{s),f{u)) c'(q; s, u) = c{q; f{s), f(u)) 

R'hiQ', s, u) = /(s), f{u)) R'^{q; s, u) = /(s), /(m)) 


(49) 

(50) 
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The striking consequence of this fact is that time here is just an arbitrary parameter. Reparametrization invariance 
plays the role of a gauge symmetry: if we want to obtain physical observables we need to fix the gauge. In this way 
we can reduce the degrees of freedom contained in the equations. A way to do it is to consider the equations in the 
regime where TTI and FDT hold. If we impose both of these conditions we get a unique relation that is given by 


0 = c(g,s) - h{q,s) + p 


h(g;s)Aco(g) + Aho(g)c(g;s) - / <lzh{q]s - z)c[q,z) + h{q-,s)c{q-,Q) 


= Wq[h]- p / dz h{q, z)[c{q, s - z) - c{q; s)] 


where 


Wg[h] = c{q; s) - h{q; s) + p [h{q; s)Aco(g) + c{q; s)Aho{q) + c{q; 0)h{q; s) - (h{q, s) - h{q, 0))c(gf, s)] . 


(51) 


(52) 


We immediately note that this equation is nothing but a MCT equation where the MCT kernel is replaced by the 
direct correlation function. This has the consequence that different modes q in the system are coupled as it should 
be since the direct correlation function can be expressed in terms of a series expansion in h{q, t) 




(- 1 )^ 


n 


d^fci 

(27r)-^ 


—d^h(fci; t) ... h{kn-i,t)h{q - ki - ... - kn-i;t) . 


(53) 


From this equation we can obtain the mode-coupling exponent parameter Amct- This quantity encodes for the 
dynamical exponents that characterize the approach and the departure of the dynamical correlation functions from 
their plateau value. The schematic way to obtain this quantity is to expand the dynamical equations using 

Ggit) = Ako{q)t’^ + SGg{t) 


hiq',t) — h{q;0) + Gq{t) 
where kyy^q) is the zero mode eigenvector [21]. In this way we get 

_r2(l-p&) 

Amct = “ 


{i+h{x)r 


(54) 


(55) 


r{l + 2b) 2pJ^k^{q){l-pAc{q))^' 

This result has been derived also from a different perspective in [21, 22]. Moreover we can also use the dynamical 
equations in the aging regime where we have done a quench of the system from a high temperature configuration 
down to a temperature lower than the MCT one. Because we are in the aging time window, we can set to zero the 
term h{q; s,0)c{q;u,0) and the dynamical equations become 


h{q] s, u) = c{q; s,u) + p 


1 /■“ 

h{q;s,u)Ac{q) + Ah{q)c{q]s,u) + ^ / dzRc{q;u, z)h{q; s, z)+ 

P Jo 


1 

+ ^ J dzRh{q;s,z)c{q-,z,u) 


(56) 


Rhiq, s, u) = Rciq; s,u) + p 


Rh{q;s,u)Ac{q) + Ah{q)Rc{q;s,u) +-: / dzRh{q; z,u)Rc{q; s, z) 


We can now consider the aging parametrization for the correlation functions 


h{q;s, 

u) 

= h{q; 

;) 

Rh{q;s, 

u) 

= -R-h 

S 

(«• 

c{q;s, 


( 

u\ 

u) 

= c[q; 

s) 

Rc{Qj 

u) 

= -Re 
s 

(?; 


and setting X = u/s, the equations become 


h{q;X) =c{q;X) + p 


h{q;X)Ac{q) + Ah{q)c{q;X) + ^ [ 

P Jo 


U; ^ A)+ 


1 

+- / dX'nhiq;X')c 
P Jo 




y/ \ sgn(A-A')' 

a] 


(57) 


(58) 


(59) 
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T^h{q,>^) = '7^c(<?; A) + p 


T^hiq; A)Ac(g) + Ah{q)TZc{q; A) + 



dA' 


T^h 



7^e(g;A') 


(60) 


By using the quasi-fluctuation dissipation ansatz 


TZhiqiX) = Px^h{q;X) 
T^c{q\X) = l5x^c{q]X) . 


(61) 


we can obtain that the value of x is fixed by the marginal stability condition according to which the dynamical 
equations must have a zero mode that in replica theory is called the replicon [3, 21, 22]. In this way all the off- 
equilibrium dynamics follow closely the one of the p-spin glass model. 


2. MCT from BPD 

It has been shown in [24] how to construct a consistent truncation scheme of the replicated BBGKY hierarchy in 
order to obtain from replicas the equation of the non ergodicity parameter that has been derived within MCT. In 
what follows we want to go beyond the non-ergodicity factor to obtain the whole MCT dynamical equation in the 
long time limit. We can do this exactly on the same lines as we did in the HNC approximation scheme. In this case 
the closure scheme is provided giving the non-diagonal elements of the replicated direct correlation function Ca^^{q) 
in terms of the static direct correlation function Co{q): 

Cai^ik) = j dqV{k,q)hap{q)hai3{k-q) (62) 

where the V{k, q) is the Mode Coupling vertex function 

V{k,q) = [k • (qco(g) -f (k- q)co(fc - g))]^ (63) 

which is independent of the replica indexes. At this point we use again the mapping of replica indices on pseudo time 
indices a —>■ (a = I, t), /? —>■ (a = I, s = 0) with t > 0. Within this scheme the direct static correlation function is 
supposed to come from equilibrium and this is the only regime we can have access to. By using TTI and FDT we get 

c{k,t) = J dq V{k,q)h{q,t)h{k — q,t). (64) 

Plugging this equation inside (44), after some simple algebra, we get the MCT equations derived by Cotze [25]. 


VI. CONCLUSIONS 

In this paper we have reviewed the construction of the Boltzmann pseudodynamics and presented some new results 
for spin glasses and liquid theory. Among them we have: 

• A close relation between the response function and the clone correlation function, which shows analytically for 
the first time how anomalous response requires non trivial clone correlations. 

• The derivation of dynamical equations for spherical models that avoid the use of the replica method. This method 
can be generalized out of mean-field to obtain a hierarchical system of equations for multibody correlation and 
response functions. 

• The results of explicit integration of the equation of motion in spherical spin-glass models, confirming the 
asymptotic analysis of the long chain limit and showing that this limit is achieved in relatively short chains. 

In addition we discussed the derivation of dynamical Ornstein-Zernike equations suggested by the formalism and we 
have showed that they have a structure that generalizes the one of the Mode Coupling equations. These equations 
can be closed using schemes borrowed from equilibrium liquid theory. We showed that if the Szamel’s closure scheme 
is applied one recovers the C otze MCT equation. An alternative is the HNC approximation which allows in principle 
a quantitative description of aging phenomena in supercooled liquids. 
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Within BPD, all the available approximations allowing to describe long time aging dynamics coherently confirm 
the original analysis of simple mean-field spin glass models, in particular effective temperatures associated to mutual 
equilibration of slow degrees of freedom naturally emerge and are interpreted. We believe that the principle of quasi- 
equilibrium configuration space exploration formalized by Boltzmann psendodynamics go beyond the approximations 
and is at the heart of a description of slow dynamics in terms of effective temperatures. 
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